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SYMBOLS 


Alphabetical  FORTRAN 

Character  Character 


AM  AM  length  of  pallet  or  load 

DE  DE  arbitrary  interval  of  time  between 

successive  points  in  solution  of  a 
differential  equation 


g 


acceleration  due  to  gravity 


G 


center  of  gravity 


HG 


HG 


height  of  center  of  gravity  of  load 
above  floor 


HH  HH  height  of  load  if  rectangular;  other¬ 

wise,  length  of  perpendicular  from 
any  point  on  upper  side  of  load  above 
floor 


I 


G 


L 


moment  of  inertia  of  load  about  hori¬ 
zontal  axis  through  center  of  gravity, 
perpendicular  to  line  of  flight 

sill  of  ramp  in  horizontal  (i.  e.  ,  up) 
position 


q 


Q 


distance  between  mid-point  of  under¬ 
side  of  load  and  sill  of  ramp,  meas¬ 
ured  positively  when  mid-point  is  rear 
of  sill 


q  QD 


first  derivative  of  q  with  respect  to 
time 


t 


time 


W 


weight  of  load  in  gravitational  units 


v 


Alphabetical 

Character 


FORTRAN 

Character 


X 


St 


X 


G 


y 


G 


e 


0 


e 


XST 


XG  YG 


T  or  TH 


TD  or  THD 


TDD 


distance  from  bulkhead  of  forward 
face  of  load  or  forward  edge  of 
pallet  when  motion  starts 

coordinates  of  center  of  gravity  of 
load  with  respect  to  horizontal  and 
vertical  axes  through  L 

x^  measured  positively  rearward 
y^  measured  positively  upward 

angle  through  which  load  has  turned 
from  original  horizontal  position; 
denoted  by  T  in  two-degree  phase 
and  by  TH  in  three-degree  phase 

first  derivative  of  0  with  respect  to 
time;  denoted  by  TD  in  two-degree 
phase  and  by  THD  in  three-degree 
phase 

second  derivative  of  6  with  respect 
to  time 


M 


YU 


coefficient  of  sliding  friction  between 
floor  and  load  or  pallet 


R  force  of  reaction  of  sill  on  load  or 

pallet  in  gravitational  units 

DIST  distance  traveled  by  load  at  constant 

acceleration  during  one-degree-of 
freedom  phase 
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SUMMARY 


A  study  is  made  of  the  extraction  of  loads  by  parachute  from  the  CV-7 
Caribou  aircraft.  Two  FORTRAN  programs  are  included  (see  Appendix 
II)  to  show  the  calculation  of  the  maximum  safe  envelope  for  the  loads; 
variouf  conditions  with  a  wide  range  of  adjustable  parameters  are  con¬ 
sidered.  In  the  programs,  the  parameters  may  be  set  as  data  to  simulate 
any  values,  such  as  those  for  the  extractive  force  of  the  ejection  para¬ 
chute,  for  the  coefficient  of  sliding  friction  between  the  floor  and  the  load, 
and  for  the  length  of  the  pallet  on  which  the  load  is  mounted. 

The  conclusions  reached  are  negative  in  character,  but  they  can  be  of 
value  as  a  basis  for  further  studies.  The  load  tips  so  little  at  the  sill  of 
the  floor  with  the  ramp  up  that  the  maximum  safe  envelope  departs  only 
slightly  from  the  rectangular  form.  It  is  believed  that  the  effects  of  gust 
disturbances  and  of  possible  jamming  in  the  rail  restraint  and  release 
system  should  be  investigated,  since,  from  time  to  time,  these  influences 
will  rause  the  load  ^^o  tip  through  angles  of  inclination  of  a  higher  order  of 
magnitude  than  those  due  to  dynamical  considerations. 
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CONCLUSIONS 


Computations  for  this  study  were  based  on  a  216-inch-long  pallet,  which 
is  the  longest  rectangle  contemplated  and  the  most  likely  to  reveal  the 
effects  of  tipping.  In  no  case  did  the  four  safe  envelopes  computed  show 
a  departure  of  more  than  7  inches  between  the  minimum  and  the  maximum 
height  pt  rmissible.  Therefore,  it  is  concluded  that: 

1.  Further  study  should  include  the  extent  of  gust  disturbance  on 
load  trajectory.  (This  study  would  be  in  addition  to  carrying 
out  experimental  work  to  determine  photographically,  or  other¬ 
wise,  the  trajectory  of  such  a  pallet.) 

2.  The  effects  of  jamming  in  the  rail  restraint  and  release  system 
should  be  investigated  by  means  of  a  nonexperimental  method: 
the  FORTRAN  data  parameters  could  be  adjusted  so  that  the  load 
would  have  a  velocity  rearward  of  almost  zero  when  tipping  starts. 
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.  INTRODUCTION 


The  study  discussed  in  this  report  is  a  continuation  of  one  presented  in 
TRECOM  Technical  Report  63-71,  "Study  of  Design  Factors  in  Air  De- 
li\t“ry  for  CV-7  Caribou  Aircraft",  which  was  published  by  the  U.  S. 

Army  Transportation  Research  Command  at  Fort  Eustis,  Virginia,  in 
December  1963.  The  earlier  report  dealt  with  the  motion  of  containers 
loaded  in  the  cargo  compartment  of  aircraft  when  the  containers  are  ex¬ 
tracted  by  parachute.  In  this  report,  consideration  is  given  to  the  effects 
of  friction  between  the  load  and  the  floor  of  the  cargo  compartment  and  to 
the  downward- sloping  configuration  of  the  roof  over  the  double-hinged 
door  at  the  rear  of  the  cargo  compartment:  also  in  this  report,  it  is  as¬ 
sumed  that  the  force  of  extraction  drops  to  zero  as  soon  as  tipping  starts 
at  the  sill  of  the  ramp. 

In  the  previous  report,  the  force  of  extraction  was  assumed  to  be  1.5 
times  that  due  to  gravity,  a  value  which  is  considerably  higher  than  that 
indicated  by  experiments  conducted  at  the  Air  Force  Flight  Test  Center, 
Edwards  Air  Force  Base,  California,  where  it  was  found  that  when  the 
parachute  was  unopened,  the  force  built  up  from  zero  value  to  a  maximum 
which  lay  between  .  3  and  .  8.  * 

These  much  lower  values  for  the  extracting  force  result  in  the  attainment 
of  increased  angles  of  inclination  of  the  load  to  the  horizontal,  and  this, 
in  turn,  produces  contours  for  maximum  safe  envelopes  which  are  no 
longer  uniform  in  height.  The  lower- the  extracting  force  and  the  greater 
the  coefficient  of  sliding  friction,  the  more  marked  this  effect  becomes. 
Four  contours  for  maximum  safe  envelopes  have  been  calculated  for  the 
216-inch-long  pallet  for  typical  values  of  such  parameters  as  the  coeffi¬ 
cient  of  friction,  the  extracting  force,  the  height  of  the  roof  from  the  sill. 


*Univer8al  Cargo  Handling  System  for  C-130  Aircraft,  Technical 
Documentary  Report  No.  62-27,  Air  Force  Flight  Test  Center,  Edwards 
Air  Force  Base,  California. 
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and  the  starting  station  of  the  pallet.  The  contours  of  three  of  the  enve¬ 
lopes  are  shown  graphically  in  Appendix  I.  The  results  for  two  of  the 
four  envelopes  differed  so  slightly  that  one  graph  represents  both  enve¬ 
lopes;  the  only  two  differences  between  the  two  envelopes  are  shown  with 
the  graph. 

This  report  also  shows  the  need  to  determine  experimentally  the  radius 
of  gyration  around  a  horizontal  axis  through  the  center  of  gravity  of  a 
load  such  as  the  105mm  howitzer  mounted  on  a  216-inch  pallet. 

In  TRECOM  Technical  Report  63-71,  the  general  equations  of  motion  of 
the  load  during  the  two-degrees-of-freedom  phase  were  derived  by 
Lagrange's  method.  Owing  to  the  presence  of  friction,  which  makes  the 
system  nonconservative,  Lagrange's  method  was  not  followed  in  the 
present  study.  Instead,  the  forces  of  reaction  and  of  friction  have  been 
included  as  external  forces  on  the  load.  One  angular  momentum  equation 
and  two  linear  acceleration-force  equations  are  derived;  these  equations 
yield  two  differential  equations  for  q  and  0  and  also  one  for  R,  the 
I  eaction. 

Appendix  IV  has  been  added  to  show  how  the  maximum  safe  height  of  a 
rectangular  load  can  be  calculated.  FORTRAN  Program  Number  Three 
has  been  included  in  this  appendix  to  show  the  calculation  for  rectangular 
loads  starting  from  any  station.  The  data  shown  at  the  foot  of  FORTRAN 
Program  Number  Three  refer  to  a  54-inch  container  starting  from  station 
291  as  in  Loading  Plan  2. 

In  the  analysis  of  the  one-,  two-,  and  three-degrees-of-freedom  motion, 
differential  equations  must  first  be  obtained  for  the  generalized  coordi¬ 
nates  that  define  the  configuration  of  the  load  at  any  time;  then,  these 
values  of  the  coordinates  are  used  to  find  whether  any  given  portion  of 
the  load  hits  the  horizontal  or  the  downward- sloping  portion  of  the  roof. 

A  general  program  is  constructed  that  can  handle  the  two-  or  three- 
degrees-of-freedom  phases  of  the  motion  and  that  during  each  phase  car 
examine  whether  either  portion  of  the  roof  is  struck.  However,  in  the 
majority  of  cases  it  is  quite  unnecessary  to  use  the  three-degrees-of- 
freedom  portion  of  the  program;  if  the  load  clears  the  two  portions  of  the 
roof  during  the  two-degree  phase,  it  clears  them  also  in  the  three-degree 
portion  of  the  motion.  For  this  reason,  FORTRAN  Program  Number  Two 
was  constructed;  it  simply  examines  the  two-degrees-of-freedom  phase 
and  gives  a  message  on  the  typewriter  as  to  whether  or  not  the  roof  is 
struck. 

In  both  FORTRAN  Programs  Numbers  One  and  Two,  the  horizontal  and 
sloping  portions  of  the  roof  are  considered.  The  height  of  the  horizontal 
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portion  of  the  roof  is  read  in  as  data  and  is  normally  considered  to  be  76 
inches  (8Z  inches  as  shown  on  the  manufacturer's  drawings  less  a  6-inch 
safety  margin).  The  equation  in  FORTRAN  Program  Number  One  for  the 
sloping  part  of  the  roof  is 

8.  lx  +  146.  3y  -  10884. 7  =  0  (1) 

in  relation  to  the  x  axis,  which  extends  horizontally  through  the  sill  of  the 
ramp,  L,  and  the  y  axis,  which  extends  vertically  through  L.  The  x  axis 
is  measured  positively  to  the  rear  and  the  y  axis  positively  upward. 

The  left-hand  member  of  this  equation  is  called  SLOP;  its  being  negative 
indicates  that  the  point  represented  by  (x,y)  is  on  the  same  side  of  the 
sloping  roof  as  the  sill,  L;  when  SLOP  is  positive,  point  (x,y)  is  on  the 
side  remote  from  L.  In  other  words,  when  SLOP  is  negative,  point  (x,y) 
has  not  touched  the  roof;  when  it  is  positive,  the  roof  has  been  struck. 

In  the  latter  case,  FORTRAN  Program  Number  Two  terminates  with  the 
message  "STRIKES  ROOF".  Otherwise,  the  message  "DOES  NOT  STRIKE 
ROOF"  appears. 

In  FORTRAN  Program  Number  One,  the  equation  for  the  sloping  portion 
of  the  roof,  which  is  indicated  in  equation  (1),  does  not  provide  a  safety 
margin  of  6  inches.  As  stated  before,  a  6-inch  safety  margin  below  the 
horizontal  portion  of  the  roof  is  provided.  Therefore,  if  a  6-inch  safety 
margin  is  desired  under  both  parts  of  the  roof,  SLOP  is  replaced  by  the 
following  equation: 

SLOP  =  8.  95xCC  +  MOyCC  -  9695. 

Some  of  the  computer  runs  have  been  made  with  this  modified  value  for 
SLOP,  and  the  results  can  be  compared  with  those  corresponding  to  equa¬ 
tion  (1).  The  same  margin  of  safety  should  be  provided  throughout.  The 
load  will  be  moving  faster  when  passing  under  the  sloping  part  of  the  roof, 
so  that  the  consequences  of  striking  the  roof  would  be  more  serious. 

In  FORTRAN  Program  Number  Two,  the  modified  value  of  SLOP  has  been 
inserted. 


ONE-DEGREE-OF-FREEDOM  MOTION 


Figure  1  shows  the  load  moving  down  the  floor  of  the  cargo  compartment 
of  the  aircraft  under  the  action  of  the  following  forces: 
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P  due  to  the  extracting  force  of  the  parachute: 

W  due  to  the  force  of  gravity  of  the  load  itself  acting  through  the 
center  of  gravity,  G,  of  the  load; 

R  due  to  the  reaction  of  the  floor  on  the  load; 

due  to  friction  of  the  floor  on  the  load. 


C 


Figure  1.  Sketch  Illustrating  One- Degree-of- 
Freedom  Motion. 

The  two  coordinates,  and  0  ,  define  the  position  and  configuration  of 

the  load  with  respect  ic  the  aircraft.  Following  is  a  theorem  of  wide 
applicability  in  the  mechanics  of  a  system  of  particles  (which  is  therefore 
applicable  to  a  rigid  body)  that  is  used  to  determine  the  horizontal  and 
vertical  components  of  the  acceleration  of  G:  the  center  of  gravity  of  a 
system  of  particles  moves  as  though  all  the  mass  and  all  the  external 
forces  were  concentrated  at  the  center  of  gravity. 

When  the  forces  are  resolved  horizontally, 

p  -  yR  s  —  (x  -  4^  cos  0  HG  sin  0), 

g  dt2  A  2 

the  quantity  in  parentheses  being  equal  to  the  x  coordinate  of  the  center 
of  gravity. 
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If  forces  are  resolved  vertically, 


W 


R  = 


w  ^  ^ 

g  dt^  g  dt2^  2 


sin  6  ♦  HG  cos  0) 


To  find  the  angular  acceleration  of  the  load  about  a  horizontal  axis 
through  G,  a  general  theorem  in  the  dynamics  of  a  rigid  body  is  again 
referenced:  the  rate  of  change  of  angular  momentum  about  any  axis 
through  the  center  of  gravity  is  equal  to  the  moment  of  the  external 
forces  about  that  axis.  Hence, 


yR(-^  sin  0  +  HG  cos  0) 


AM 

R(-i^  cos  0  -  HG  sin  0) 


-  PC^  sin 


0  (HG  -  HD)  cos  0]  = 


g 


When  the  differentiations  are  performed,  these  equations  become 

P  -  uR  =  -Cx  t  0^(Ar  cos  0  -  HG  sin  0)  +  0(4Ji  sin  0  ♦  HG  cos  0)1  (2) 
g  A  2  2 


W  -  R  =  -[02(M  sin  0  HG  cos  0)  ♦  0(HG  sin  0  -  cos  0)]  (3) 

g  2  2 


uR(^  sin  0  HG  cos  0)  -  R(^  cos  0  -  HG  sin  0) 

-  P[^  sin  0  +  (HG  -  HD)  cos  0]  *  0i^.(4) 

^  e 

In  order  to  find  the  conditions  under  which  the  load  will  tip  about  the 

corner  A,  let  0=0  and  0=0.  Also  set  P  =  W  PRO  and  I  -  *  4HG2 

G  12 

by  the  well-known  formula  for  the  moment  of  inertia  of  a  rectangular 
parallelepiped  about  an  axis  through  its  center  of  gravity  parallel  to  an 
axis  of  figure.  If  these  substitutions  are  made,  the  equations  become 
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W  PRO 


(5) 


aUG-  ♦  Ru  +  X.  -  = 
R  g 


.0^  -  +  R  =  W 
2  P, 


^  _  yHQ)  _  .y  pro(HG  -  HD)  . 

g  12  2 


Solving  by  Cramer's  rule  for  0  , 


e  s  g 


PRO 

1 

.PRO(HG  ■ 


^  am 

HD):  ^  -  uHG 


HG 

2 

(AM2  t  4HG2) 

12 


u 

1 


AM 

'Y  " 


w/g 

0 

0 


’•^/g 

0 


(6) 

(7) 


Tiie  denominator  is  negative  for  all  values  of  AM,  PRO,  HG.  and  HD  that 
will  encountered  in  practice.  Hence,  the  initial  value  of  0  is  positive 
if  the  numerator  is  negative,  that  is,  if 


uHG  ♦  PRO  HD  >  -^  +  PROHG. 


(8) 


Provided  this  condition  is  not  satisfied,  0  and  its  first  two  derivatives 
remain  zero  until  tipping  starts  at  the  rear  end  of  the  loading  ramp  during 
the  two-degrees-of-f reedom  motion.  As  long  as  0  =  0,  equations  (4), 

(5),  and  (6)  reduce  to 


and 


Ru  +  X.  -  s  W  PRO 
A  g 

R  =  W, 


which  yields  x  =  g(PR0  -  u)  otherwise  evident.  (9) 

This  value  of  x/^  is  used  in  line  16  of  FORTRAN  Program  Number  One 
and  in  line  15  of  FORTRAN  Program  Number  Two  to  give  the  value  of 
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QD(1),  that  is,  the  initial  linear  velocity  rearward  of  the  load  at  the  start 
of  the  two-degrees-of-freedom  motion.  The  load  will  have  enjoyed  this 
constant  acceleration  since  starting  from  rest  with  corner  B  at  distance 
Xgj  from  the  forward  bulkhead  until  tipping  starts  when  B  is  at  distance 

D1  -  [.^  -  Q(l)]  from  the  bulkhead.  The  difference,  which  is  called 

DIST  in  the  FORTRAN  programs,  between  these  two  distances  gives  the 
length  of  the  path  traveled  at  constant  acceleration  g(PRO  -p). 


TWO-DEGREES-OF-FREEDOM  MOTION 


As  stated  earlier,  it  is  assumed  in  this  study  that  the  force  of  extraction 
of  the  parachute  becomes  zero  the  moment  that  tipping  starts  at  the  rear 
end  of  the  loading  ramp.  Figure  2  shows  how  the  configuration  of  the 
system  is  defined  by  the  two  generalized  coordinates,  q  and  0  .  The 
three  external  forces  acting  on  the  load  are  also  shown:  the  force  of 
gravity  vertically  downward  through  G;  the  reaction  of  the  sill  on  the 
underside  of  the  load,  R,  acting  normally  to  the  underside;  and  the  force 
of  friction,  piR,  tangentially  to  the  underside. 


c 


Figure  2.  Sketch  Illustrating  Two- Degrees-of 
Freedom  Motion. 
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If  the  forces  are  resolved  horizontally. 


-uR  cos  6  t  R  sin  0  =  0  ♦  HG  sin  0)  “  0  -  q0 

•*  •  <>  •• 

sin  0  -  qe  sin  6  -  q0^  cos  0  -  q®  sin  0  -  HG0^  sin  0  ♦  HG0  cos  0).{1O) 
If  the  forces  are  resolved  vertically, 

d2 

W  -  R  cos  0  -  uR  sin  0  s  --^(HG  cos  0  -  q  sin  0) 

W  •  o  ••  ** 

=-|(-HG0^  cos  0  -  HG0  sin  0  -  q  sin  0 

-  2q6  cos  0  +  q0^  sin  0  -  q0  cos  0),  (11) 


since  vertical  coordinates  are  measured  positively  upward. 


Taking  moments  about  the  center  of  gravity,  G, 


WRHG  +  Rq  ♦  Iq  H  , 


(12) 


Set 


I 


G 


,„AM2  ♦  4HG2, 
W( - ) 

12 


and  arrange  equations  (10),  (1  1),  and  (12)  as  equations  in  the  second 
derivatives  of  q  and  6  ;  then  in  R, 


#• 

q  cos  0  f  0(HG  cos  0  -  q  sin  0) 


-  2q0  sin  0  +  0^(q  cos  0  +  HG  sin  0),  (13) 


10 


and 


q  sin  0  ♦  e(HG  sin  0  ♦  q  cos  0)  t 


_ R _ 

Wg  (cos  0  f  u  sin  0) 


=  g  -  e2(HG  cos  0  -  q  sin  0)  -  2q0  cos  0,  (14) 


and 


q 


„  •vAM2  ♦  4HG2 
0  +  0( - rs - 


Rg(q  ♦  P  HG)  _  ^ 
W  ~  ^ 


(15) 


From  this  point  on,  both  in  the  analysis  and  in  the  FORTRAN  programs, 
R  is  replaced  by  R. 

W 

Equations  (13),  (14),  and  (15)  take  the  following  simplified  forms: 

qAl  +  0A2  ♦  RA3  =  A4  (16) 

qBl  f  0B2  t  RB3  3  B4  (17) 

qO  ♦  ’k2  t  RC3  =  0  .  (18) 


The  following  substitutions  can  be  made: 


A1  s  cos  0,  A2  s  HG  cos  d  -  q  sin  0,  A3  -  g(u  cos  0  -  sin  0),  and 
A4  s  2q0  sin  0  ♦  0^(q  cos  0  ♦  HG  sin  0). 


B1  s  sin  0,  B2  s  HG  sin  0  +  q  cos  0,  B3  =  g(cos  0  u  sin  0),  and 
B4  s  g  -  02(hg  cos  0  -  q  sin  0)  -  2  q0  cos  0. 


C2 


AM^  »  4  HG^ 
12 


and  C3  =  g(-q  -  uHG)  . 


Then  equations  (16),  (17),  and  (18)  can  be  solved  for  q,  6  ,  and  R  by 
Cramer's  rule.  The  10  quantities  A1  to  C3  retain  their  identities  in  the 
FORTRAN  programs.  The  values  obtained  for  R  must  always  be  positive; 
this  has  been  tested  in  all  instances. 


1 1 


In  FORTRAN  Program  Number  One  (see  Appendix  II),  immediately  pre¬ 
ceding  the  "51  Continue"  statement  that  closes  the  first  DO  loop,  the 
following  two  "IF"  statements  appear: 

IF  (Q(K)-AM/2.  )  12,  13,  45 
12  IF  (Q(K  +  l)-AM/2.  )  51,  14,46  . 

These  statements  direct  the  program  to  four  alternatives,  which  in  each 
case  calculate  by  simple  interpolation  the  values  of  q,  0,  and  0  ;  these 

AM 

values  correspond  to  q  =  — >  that  is,  to  the  instant  when  corner  B  is  at 

L,  the  sill  of  the  ramp.  These  final  values  of  q  and  of  0  and  of  their 
first  derivatives  for  the  two-degrees-of-freedom  motion  are  used  to 
obtain  initial  values  of  coordinates  and  their  velocities  in  the  three- 
degrees-of-freedom  phase. 


THREE-DEGREES-OF-FREEDOM  MOTION 


Immediately  after  corner  B  of  the  load  reaches  the  sill,  L,  the  load 
assumes  three  degrees  of  freedom,  the  remaining;  three  of  the  six  free¬ 
doms  enjoyed  by  a  rigid  body  in  space  being  ankylosed,  as  explained  in 
TRECOM  Technical  Report  63-71  on  page  1.  The  only  external  force 
acting  on  the  load  during  the  three-degree  phase  is  that  of  gravity;  con¬ 
sequently,  the  center  of  gravity,  G,  executes  a  parabolic  path,  and  the 
load  continues  to  turn  with  uniform  angular  velocity  about  a  horizontal 
axis  through  G. 

Let  Xq  and  y^  be  the  coordinates  of  G  with  respect  to  horizontal  and 

vertical  axes  through  L,  x  being  measured  positively  aft  and  y  positively 

••  •• 

upward.  Thus,  x^  =  0  and  y^  =  -g,  and  on  integrating  twice  with  respect 
to  time. 


■  ‘’‘6*0 


*  ""g’o 


and 


''g 


‘''g'o  ♦  "''g>o 


-  8^ 
1" 


where  t  seconds  have  elapsed  since  corner  B  was  at  L.  Also,  in  an 
obvious  notation,  (.  .  )q  is  used  to  denote  the  value  of  a  variable  at  the 
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•  • 

instant  when  t  =  0.  The  final  values  of  q,  q,  0,  and  0  in  the  two-degrees- 
of-freedom  phase  are  denoted  by  QL,  QLD,  TL,  and  TLD. 

At  any  instant  in  the  two-degree  phase, 

=  q  cos  0  HG  sin  0 

G 

and 

•  •  •  t> 

Xg  s  q  cos  0  -  q0  sin  0  HG0  cos  0, 

Hence, 

(x  )  =  QL  cos  TL  ♦  HG  sin  TL 

G  o 

and 

(x  )  =  QLD  cos  TL  -  QL  sin  TL  TLD  +  HG  TLp  cos  TL  . 

G  o 


(19) 

(20) 


Similarly, 


Hence, 

(y  )  =  HG  cos  TL  -  QL  sin  TL,  (21) 

G  o 

and 

•  •  •  • 

y  s  -HG0  sin  0  -  q  sin  0  -  q0  cos  0, 


y  =  HG  cos  0  -  q  sin  0. 
^G 


which  in  turn  yields 

(y  )  s  -HG  sin  TL  TLD  -  QLD  sin  TL  -  QL  TLD  cos  TL,  (22) 

In  FORTRAN  Program  Number  One,  the  quantities  on  the  left-hand  sides 
of  equations  (19)  through  (22)  are  referred  to  as  XGO,  XGDO,  YGO,  and 
YGDO. 

In  TRECOM  Technical  Report  63-71,  pages  20  through  23,  the  differential 
equation 

•  • 

#  s  -SS  sin  ^ 

was  solved  to  yield  the  angle  of  inclination  of  the  load  to  the  horizontal 
during  the  three-degrees-of-freedom  phase.  However,  in  our  present 
problem,  as  previously  stated,  the  extracting  force  of  the  parachute  be¬ 
comes  zero  when  tipping  at  L  commences.  Hence,  the  load  turns  with 
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uniform  angular  velocity  around  G  during  the  three-degrees-of-fr eedcin 
motion;  and  the  angle  of  inclination  of  the  load  to  the  horizontal,  which 
is  denoted  by  TH  in  FORTRAN  Program  Number  One,  is  found  at  succes¬ 
sive  intervals  of  time  by  the  following  FORTRAN  equations: 

TH(J  +  1)  =  TH(J)  +  DE*THD(J) 

THD(J  +  1)  =  THD(J). 

In  the  FORTRAN  programs  show’n  in  TRECOM  Technical  Report  63-71, 
the  following  expression  was  incorporated  for  the  moment  of  inertia  of 
the  load  around  a  horizontal  axis  through  its  center  of  gravity: 

AM2  -f  4HG2 
*  12 

Using  this  equation  is  tantamount  to  assuming  that  the  load  is  a  rectan¬ 
gular  parallelepiped  of  length  AM  and  height  2HG.  While  this  is  a  fair 
approximation  for  cargo  loaded  into  a  container,  the  expression  is  far 
from  exact  for  an  irregular  object  such  as  a  105mm  howitzer,  which  has 
much  of  its  weight  concentrated  around  its  center  of  eravitv.  For  this 
reason,  in  FORTRAN  Programs  Numbers  One  and  Two  the  quantity  C2, 
which  is  the  moment  of  inertia  in  equations  (16)  through  (18),  is  reduced 
by  an  arbitrary  proportion  that  in  the  sample  program  in  Appendix  II  is 
taken  as  0.6.  While  this  is  an  arbitrary  and  inexact  method  of  estimating 
the  moment  of  inertia,  it  is  the  best  way  available  until  actual  experiments 
are  conducted  to  determine  this  quantity.  The  moment  of  inertia  of  such 
a  load  could  be  determined  by  timing  its  periods  of  oscillation  around  a 
ser  ies  of  axes  that  are  parallel  to,  and  at  varying  distances  from,  the 
load's  center  of  gravity. 

The  quantity  C2  has  been  further  modified  to  correspond  to  a  proportion 
of  .4  times  that  for  a  uniform  rectangular  load  in  certain  of  the  computer 
results;  in  this  case  the  angle  of  inclination,  G  ,  of  the  load  increases 
more  rapidly.  In  all  of  the  graphical  results  of  the  computer  runs  and  in 
the  corresponding  tables  of  the  values  of  HH  and  AN,  the  proportion 
assumed  for  modifying  C2  is  clearly  stated. 
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APPENDIX  I 


MAXIMUM  SAFE  ENVELOPES 


DI3TANCL  rkOM  fORWARD  HLI. 

;i!T  ABOVC  riooR 

PARAMLTLRS 

rAiX  CT  b.iAh  (In.) 

(  In.  ) 

Acce ler.1T  ive  force  is  ,4  tines  tnat  of  pravity. 

J 

7j 

1  Coefficient 

of  friction  is  .03, 

IH 

73 

distance  f  ron  buikrie.id  when  load  starts  novinp  with 

Jb 

74 

forward  face  is  7d  inrhes. 

74 

Lenptfj  of  pallet  on  which  load  is  mounted  is  ?1^ 

17 

73 

1  inches. 

‘JO 

73 

Center  of  gravity  of  load  is  ?J  inches  af>ove  nid- 

108 

17 

point  of  pallet. 

i2t) 

17 

safety  nvargin  provided  l^elow  horizontal  part  ol  roof 

144 

71 

is  b  inches;  no  safety  margin  is  provided  below 

162 

70 

sloping  part  of  roof. 

180 

69 

Moment  of  inertia  ol  loa<I  is  ,b  times  that  of  rec- 

198 

68 

tanpular  loatl  of  same  dimensions. 

?lfa 

67 

90- 

•  €) 

1 

^  1 

i  c 

MaxiinurTi  height  of  load  at  45  inches  from  leading  edge  74  inches 

0  :: 

JJ  u 

< 

Minimum  height  of  load  at  rear  edge  67 

inches 

X  E  30- 

- 

QC  O 

1  i  1  1 

4  1 

1 

< 

1  f  1  f 

)  30  60  90  120 

1  t 

150  160 

210 

Distaiue  K  r  >m  Furwarc  Kdg* 

of  Load  (in.  ) 
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DISTANCL  FFOH  PORi^ARO 
FACE  or  LOAD  (In.) 


HEIGHT  AitOVE  FLOOR 
(In.) 


Accel«rativ«  fore*  is  .2^  tiass  that  of  gravity. 

Coefficient  of  friction  is  .Ob. 

Distance  froa  bulkhead  when  load  starts  aoving  with 
f'^rward  face  is  78  inches. 

of  pallet  on  which  load  is  aounted  is  216 
inches. 

Center  of  gravity  of  load  is  20  inches  above  aid- 
point  of  pallet. 

Safety  aargin  provided  below  horisontal  part  of  roof 
is  b  inches;  no  safety  aargin  is  provided  below 
sloping  part  of  roof. 

Moaent  of  inertia  of  load  is  .6  tiaes  that  of  rec¬ 
tangular  load  of  saae  diaensions. 


PARAMETERS* 


Distance  From  Forv^ard  Edge  of  Load  (in.  ) 


DISTAHCE  FROM  FORSAKO 
FACE  or  LOAD  (In.) 


HEIGHT  ABOVE  FLOOR 
(In.) 


Accelerative  force  is  .25  tiaes  that  of  gravity.  0  71 

Coefficient  of  friction  is  .05.  18  73 

Distance  froa  bulkhead  when  load  starts  aoving  with  3b  74 

forward  face  is  22  inches.  54  74 

Length  of  pallet  on  which  load  is  aounted  is  216  72  74 

inches.  90  73 

Center  of  gravity  of  load  is  20  inches  above  aid-  108  73 

point  of  pallet.  1?6  72 

Safety  Mrgin  provided  below  horizontal  part  of  roof  144  7X 

is  6  inches;  no  safety  aargin  is  provided  below  162  70 

sloping  part  of  roof.  180  69 

Hoawnt  of  inertia  of  load  is  .6  times  that  of  rec-  198  68 

tangular  load  of  saae  diaensions.  216  67 

*  Results  froa  the  fourth  envelope  differed  in  only  the  following  respects:  (1)  the  moaent  of 
inertia  for  the  fourth  envelope  was  .4  (rather  than  .6)  times  that  of  a  rectangular  load  of 
the  saae  diaensions;  an ^  (2)  when  there  was  a  distance  of  36  inches  from  the  forward  face  of 
the  load,  the  height  above  the  floor  was  73  (rather  than  74)  inches. 


Disf  iiu f  From  Forward  Edge  of  Load  (in.  ) 
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APPENDIX  II 


FORTRAN  PROGRAMS 


PROGRAM  NUMBER  ONE 

C  CARIBOU  WITH  FRICTION 

DIMENSION  T( 100) ,TD( 100) tTDD(  100 ) . 0 ( 100 ) .QD ( 1 00 ) »ODD( 100)  #R( 100) 
DI'^FNSION  TH(45)»THD(45)  »THDD(45) 

READ  ItPRO 

READ  l.YU 

READ  1*DE 

READ  2#D1 

READ2»XST 

RFAD2#AM 

READ  2*HH 

READ  3#AN 

read  2»HG 

READ  2tROOF 

Q(  1)— YU*HG 

DIST«Dl-XST-AM/2.+Q( 1 ) 

QD( 1 ) -SORT! 773.0* (PRO-YU )*D I  ST) 

T(1)«.C 

TD(1)».0 

C2-(AM**2+4,*HG**2)/12. 

C2*.6*C2 
DO  51  K-1.65 
Al»COS(T(X) ) 

A2-HG*COS(T(IC)  )-(Q(K)  )*SIN(T(K)  ) 

A3»386.4*(YU*COS(T(K) )-SlN(T()0  ) ) 

A5»2.*(OD(K  )  )♦(  T0(K)  )*SIN(  T(IC»  ) 

A6»(  (TD(K) )**2)*( (Q(K) )*COS(T(K) )+HG*SIN(T(K ) ) ) 

A4« A5+A6 
B1«SIN(T(K) ) 

B2-HG*SIN(TtK)  )*(Q(K) )*(COSIT(K) ) ) 

B3«ICOS(T(K)  )+YU#SrN(T(K) )) *386.4 

B5«306.4-(  (  TD(K)  )  **2  )  *  (  HG*COS  (  T  ( K  )  )-(0(IC))*SIN(T(X.))) 
B6--2.*(0D(K;)  )*(TD(K)  )*C0S(T(IC)  ) 

B4-3A+B6 

C3»-386.4*(0(X)+YU*HG) 

DET=A1*(B2*C3-B3*C2)-A2*B1*C3+A3*B1*C2 

0DD(K)«IA4*IB2*C3-C2*B3)-A2*B4*C3+A3*B4*C2)/DET 

TDD(K)»(A1*B4*C3)/DET 

R(X)-(-Al*B4*C2+A4*Bl*C2)/DET 

XCC«HH*SIN(T(K) )-(AN-Q(K) )*COS(T(K)  ) 

YCC»HH*COS( T(K)  )♦( AN-Q ( K) ) •SIN ( T ( X)  ) 

SLOP«XCC*8.1+YCC*146. 3-10004.7 
IF  (ROOF-YCC)  61.61.53 

53  TF  (SLOO)  54.61.61 

54  CONTINUE 

0(X*1)»0(K)+DE*0D(K)*( (0E**2)/2. )*QDO(K) 

T(X+1  )»TIK)-*-DE*TD(IC)  +  (  (DE**2)  /2.  )*TDD(K) 

TD()C+1  )-TD(X)*0E*T00(K) 

QD  ( (C*  1 ) -OD  (  K )  ♦0E*O00  ( X  ) 

IF  (0(K)-AM/2.)  12.13.45 

12  IF  (0(K*1 )-AM/2.)  51.14.46 
51  CONTINUE 

13  CONTINUE 
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C  0(0  COINCIDES  WITH  SILL  OF  FLOOR 

99  0L=  AM/?, 

0LD«00(O 
TL=T ( < ) 

TLD  =  T0(O 
GO  TO  A9 
A5  CONTINUE 

C  SILL  IS  BETWEEN  O(K-l)  AND  Q(K) 

C  WE  FIND  TM  THE  TIME  FROM  INSTANT  K-1  TILL  LOAD  CLEARS  FLOOR 

98  TM  =DE*(  AM/2, -QdC-l  )  )/(Q(<)-Q(K-l  )  ) 

0L=AM/2 , 

0LD  =  QD(K-1  )+TM»ODD(»C-l  ) 

TL*  T(K-1  )  +  TM»TD (<-!)  +  ( TM»*2)*(TDD( K-1 ) )/2, 

TLD*TD(K-1  )+TM»TDD(K-l ) 

GO  TO  49 
14  CONTINUE 

C  0(K+1)  COINCIDES  WITH  SiLL  OF  FLOOR 
97  0L*a^/2, 

OLD^QD ( K+1 ) 

TL  «T(K+1) 

TLD  »TD(K+1) 

GO  TO  49 
46  CONTINUE 

C  SILL  IS  BETWEEN  OIK)  AND  0(K4-1) 

C  WE  FIND  TM,  THE  TIME  FROM  INSTANT  K  TILL  LOAD  CLEARS  FLOOR 

96  TM  =  DE*(AM/2,-0(K) )/(0(K+l)-0(K) ) 

0L=  AM/2. 

0LD«0D(K)+TM»0DD(K) 

TL=  T(,0+  TM#TO(K)  +  (  TM**2  )  *  (  TDD  (  K  )  ) /2  . 

TLD  =  TD<<)  +  TM*TDD(0 
49  CONTINUE 
TH( 1 ) »TL 
THD( 1 )«TLD 
59  DO  61  J=l,40 
TP«J 

TT’  =  DE*TP 

XGO=OL*COS( TL)+HG*SIN(  TL ) 

XGDO=QLD*COS(TL)+TLD*(-QL*SIN( TL)+HG*C0S(TL) ) 
YGO«-OL*SIN(TL)+HG*COS(TL ) 

YGDO=-OLD*SINITL )-TLD«(QL»COS( TL)+HG»SIN( TL) > 
XGa/GO+XGDO*TIM 

YG«YGC+YGDO*TIM-( ( 386 , 4* ( T I M ) **2 ) / 2 . ) 

XC»XG-AN*COS(TH( J) )+(HH-HG)*SlN(TH( J) ) 

YC»YG  +  AN*SIN(TH( J) )  +  (  HH-HG ) *COS ( TH  (  J )  ) 

PRINT  7,J,XC»YC 
PAUSE 

TH(  J+1 )  =  TH( J)+DE*THD( J) 

THD( J+1 ) *THD( J) 

61  CONTINUE 
STOP 

1  FORMAT  (F3,2) 

2  FORMAT  (F5.1) 

3  FORMAT  (F6.1) 

7  FORMAT  (14, 3X  E14.8,3X  E14.8) 

END 
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PROGRAM  NUMBER  TWO 


C 


53 

54 


51 

71 


61 

41 

1 

2 

3 

1  1 
12 


CARIBOU  WITH  FRICTION 

DIMENSION  T( 100) »TO( 100) tTOO( 100 ) *0 ( lOO) »QD ( 100 )  .ODDI  100) tR ( 100) 

READ  l.PRO 

READ  1*YU 

READ  1*DE 

READ  2#D1 

READ  2»XST 

READ  2»AM 

READ  2*HH 

READ  3. AN 

READ  2*HG 

READ  2. ROOF 

0( 1 )«-YU«HG 

DIST«Dl-XST-AM/2.+Q( 1 ) 

0D( 1 ) -SORT! 773.0* (PRO-YU )*D I  ST) 

T  (  1 )  *  •  0 
TD( 1 ) -.O 

C2«(AM«»2+4,»HG**2) /12. 

C2-.6*C2 
DO  51  K«lt99 
A1«C0S(TrK) ) 

A2»HG*C0S(T(K) )-(0(K))*SlNCT(K)) 

A3»306.4*(YU*COS(T(K) )-SlN(T(K))) 

A5=2.*(0D(K) )*(TD(K) )*SIN(T(K) ) 

A6»(  (  TD(K)  )»*2)*(  (Q(<)  )*COS(T(K)  )+HG*SIN(  T  ( 1C)  )  ) 

A4s A5*A6 
B1-S1N(T(K)  ) 

B2«HG*SIN(T{IC)  )  +  (0(IC)  )*(COS(T(K)  )  ) 

B3*(C0SIT(K)  )+YU*SINi  T(IC)  )  )  *386.4 

B5«386.4-{  (TD()C)  )  »*2  )  *  (HG*COS  (  T  (K  )  )  - 1 0  (K )  )*SIN(T(K)  )  ) 
B6=-2.*(QD(K) )*(TD(K) )*COS(T(<) ) 

B4*B5+B6 

C3*-306.4*(O(K)+YU*HG) 

DET"A1*(B2*C3-B3*C2 )-A2*Bl*C3+A3*Bl*C2 
ODD()C)«(A4*(B2*C3-C2*03)-A2*B4*C3+A3*B4*C2  )  /DET 
TDD(IC)  =  (A1*B4*C3)/DET 
XCC*HH*S1N(  T(K)  )-(AN-0(IC)  )*COS(T(K)  ) 

YCC*HH*COS(T(K) )+(AN-Q(K) )*SIN(T(X) ) 

SLOP=XCC»y.95  ♦YCC*140. 0-9695.0 
IF  (ROOF-YCC)  61»61.53 
IF  (SLOP)  54.61.61 
CONTINUE 

Q(IC  +  1  )»Q(K)-*-OE*00()C)  +  (  (DE**2) /2.  )*000(IC) 

T(K  +  1  )'T(K)+DE*TD(K)  +  ( (DE**2) /2. )*TDDIK) 

TD(X  +  1  )«TD(K)+DE*TDD(K) 

0D(K  +  1  )=0D(X)+0E*0DD(K) 

IF  (0(K)-4M/2.)  51»51t71 

CONTINUE  , 

PRINT  11 
GO  TO  41 
PRINT  12 
CONTINUE 
STOP 

FORMAT  (F3.2) 

FORMAT  (F5.1) 

FORMAT  (F6.1) 

FORMAT  (21H  DOES  NOT  STRIKE  ROOF) 

FORMAT  (13H  STRIKES  ROOF) 

END 
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APPENDIX  IV 


MAXIMUM  SAFE  HEIGHT  OF  RECTANGULAR  LOAD 


The  maximum  height  for  a  54-inch-long  container,  starting  from  rest 
with  the  forward  face  at  station  291,  as  in  Loading  Plan  2,  can  be  found 
by  using  FORTRAN  Programs  Numbers  One  and  Two,  taking  the  least 
height  corresponding  to  any  value  of  AN. 

Exactly  the  same  result  is  found  if  FORTRAN  Program  Number  Three 
(shown  below)  is  used.  First,  the  program  is  run  in  abridged  form  from 
line  1  to  the  statement  "51  CONTINUE",  with  a  print-out  for  Q(K),  which 

reveals  that  Q(16)  is  the  first  Q  to  exceed  AM,  Accordingly,  the  second 

2 

part  of  the  orogram,  that  is,  from  "51  CONTINUE"  to  "END",  is  run 
with  the  upper  value  16  for  index  K,  as  follows: 

DO  54  K  =  1,  16. 

The  message  "STRIKES  ROOF"  is  given  with  M  =  7  and  HH(M)  =  68. 

This  shows  that  67  inches  is  the  maximum  safe  height  for  a  54-inch-long 
container  starting  from  station  291. 


FORTRAN  PROGRAM  NUMBER  THREE 

FORTRAN  PROGRAM  NUMBER  THREE 


C  CARIBOU  FOR  54  INCH  LOAD 

DIMENSION  T(50) *T0(50) «TDD(50)»0(50)>OD(50)»OOD(50) 
DIMENSION  HH(30».AN(13) 

READ  1»PR0 

READ  ItYU 

RFAO  1»DE 

READ  2f01 

READ  2»XST 

READ  2»AM 

RFAO  2»HG 

READ  2»R00F 

Q(  1 1—YU»HG 

DIST-ni-XST-AM/2.-H3(  1.) 

QD(  1  )-S0RT(773.8*(PR0‘-YU)*DIST) 

T( 1 )-.0 
TD(1)-.0 

C2"(AH**2+4.»HG»*2)/12, 

DO  51  K«l|40 
A1-C0S(T(K) ) 

A2«HG*C0S(T(K)  )-(Q(li:))*SIMT(K)) 
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A3»306.4*(YU*COS(T(IC)  )-SIN(T(K)  )  » 

A5-2.*{QD(K)  )*(TO(IO)*SIN(T(K)) 

A6»(  (  TD(IC)  )**2  )*(  (U(K)  )»COS(T(K)  ) +80*5  INC  T(  K  )  )  J 

A4>A5-»>A6 

B1«SIN(T(K)  ) 

B2«HG*SIN(T(K.)  )  +  {0(K)  )*C0S(T(K)  I 
B3«{COS(T(K)  )+YU*SIN(T(K» ) I •386.4 

B5-386.4-((TD{K)  J**2)*(HG*C0S(T(K)  ) -( Q  ( K)  )*SIN{  T  (K  )  ) ) 
B6--2.*COD(<) )*(TD(K)  )*C0S(T(K) ) 

B4-B5+86 

C3«-386.4*{Q{ K)+YU*HG) 

DET-A1*(B2*C3-B3*C2)-A2*B1*C3+A3*B1*C2 

ODD  ( 1C)  •CA4*(B2*C3-C2*B3)-A2*B4*C3 +  4  3*64*02  )/DET 

TDD(IC)-(A1*B4*C3I/DET 

0(  K+1  )«0(IC)+DE*00(K  )  +  (  (DE**2)  /2.  )  *QDD(  K) 

T(K  +  l)»T(K)+DE*TD(tC)  +  (  ( DE**2 ) /2.  )  *  T  DD  ( K  ) 

TD(K  +  1  )»TD(IC)+DE*TDDC<) 

QD(K  +  1)-0D(IC)+0E*QDD(<) 

IF  (0(K)-AM/2.  )  51 .51»58 
51  CONTINUE 
58  CONTINUE 
HH( 1 )»6t.O 
DO  56  M«l»25 
ANtl)-27.0 
DO  55  J-1.10 
DO  54  IC-1.16 

XCC»HH(M)*SIN(  T(K)  ) -(  AN  (  J  ) -Q(  1C )  )*COS(T(K)  ) 
YCC»HH(M)*COS  (  T(IC))  +  (AN{J)-0(IC))*SIN(TCIC)) 
SL0P«XCC*8.95+YCC*1 40.0-9695.0 
IF  (ROOF-YCC)  61. 61.53 

53  IF  (SLOP)  54.61.61 

54  CONTINUE 

AN( J+1 )«AN( J)-6.0 

55  CONTINUE  - 

HH(M+1 )»HH(M)+1,0 

56  CONTINUE 
PRINT  11 
GO  TO  41 

61  PRINT  12 
41  CONTINUE 

PRINT  5.M.HH(M) 

STOP 

1  FORMAT  (F3.2) 

2  FORMAT  (F5.n 

5  FORMAT  ( 15. 5X  E14.8  ) 

11  FORMAT  (21H  DOES  NOT  STRIKE  ROOF) 

12  FORMAT  (13H  STRIKES  ROOF) 

END 


.40 

.05 

.01 

450.0 

291.0 

054,0 

000,0 

076.0 
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DISTRIBUTION 


U.  S.  Army  Materiel  Command 
U.  S.  Army  Mobility  Command 
U.  S.  Army  Aviation  Materiel  Command 
U.  S.  Army  Transportation  Research  Command 
U.  S.  Army  Engineer  Research  and  Development 
Laboratories 

Army  Research  Office-Durham 
U.  S.  Army  Test  and  Evaluation  Command 
U.  S.  Army  Combat  Developments  Command 
U.  S.  Army  Combat  Developments  Command 
Aviation  Agency 

U.  S.  Army  Combat  Developments  Command 
Transportation  Agency 

U.  S.  Army  Command  and  General  Staff  College 
U.  S.  Army  Transpo  tation  School 
U.  S.  Army  Aviation  School 
Deputy  Chief  of  Staff  for  Logistics,  D/A 
U.  S.  Army  Transportation  Center  and  Fort  Eusti 
U.  S.  Army  Aviation  Test  Board 
U.  S.  Army  Airborne,  Electronics  and  Special 
Warfare  Board 

U.  S.  Army  Aviation  Test  Activity,  Edwards  AFB 
U.  S.  Army  Transportation  Engineering  Agency 
NASA  Representative,  Scientific  and  Technical 
Information  Facility 

U.  S.  Army  Standardization  Group,  Canada 
Canadian  Liaison  Officer, 

U.  S.  Army  Transportation  School 
Defense  Documentation  Center 
U.  S.  Army  Research  and  Engineering  Center 
Wright-Patterson  AFB,  USA  Liaison  Detachment 
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